Enhancement of shot noise due to the fluctuation of Coulomb interaction 
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We have developed a theoretical formalism to investigate the contribution of fluctuation of Coulomb interac- 
tion to the shot noise based on Keldysh non-equilibrium Green's function method. We have applied our theory 
to study the behavior of dc shot noise of atomic junctions using the method of nonequilibrium Green's function 
combined with the density functional theory (NEGF-DFT). In particular, for atomic carbon wire consisting 4 
carbon atoms in contact with two Al(lOO) electrodes, first principles calculation within NEGF-DFT formalism 
shows a negative differential resistance (NDR) region in I-V curve at finite bias due to the effective band bottom 
of the Al lead. We have calculated the shot noise spectrum using the conventional gauge invariant transport 
theory with Coulomb interaction considered explicitly on the Hartree level along with exchange and correlation 
effect. Although the Fano factor is enhanced from 0.6 to 0.8 in the NDR region, the expected super-Poissonian 
behavior in the NDR region is not observed. When the fluctuation of Coulomb interaction is included in the 
shot noise, our numerical results show that the Fano factor is greater than one in the NDR region indicating a 
super-Poissonian behavior. 
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I. INTRODUCTION 

Quantum effects have become remarkably significant in 
nanoscale semiconductors and the traditional Boltzmann 
equation is no longer sufficient to describe transport phenom- 
ena. As proposed by Schottky in his seminar work^i parti- 
tion noise, or shot noise, is resulted from the quantization 
of charge. Accordingly, when electrons are uncorrelated, the 
classical value of Fano factor, which describes the magnitude 
of the electric fluctuation, should be one. When a Fano fac- 
tor deviates from one, it shows a signature of interactions be- 
tween current flow in different probes. It is known that shot 
noise is influenced by two crucial factors, namely, Pauli prin- 
ciple and Coulomb interaction, which coexist in electronic 
systems. Specifically, Pauli interaction can only suppress the 
Fano factor below one which corresponds to a sub-Poissonian 
case and has been confirmed convincingly by experiments i^i^ 
The Coulomb interaction, however, could either reduce the 
shot noise, or enhance it so that the Fano factor shows a super- 
Poissonian value, depending on details of mesoscopic devices. 
Hence, quantum enhancement of shot noise from the classical 
value has been the subject of growing interest in recent years 
and is explored intensivelyi^ 

In mesoscopic systems shot noise is very important since 
it provides abundant information about transport properties 
of conductors, such as kinetics of electronsjiS distributions 
of energy," and correlations of electronic wave functions. 
In addition, people observed experimentally a shot noise en- 
hancement in the negative differential resistance (NDR) re- 
gion with or without the magnetic field in the tunneling 
structure Various mechanisms have been proposed to gen- 
erate a NDR, including enhancement of tunnel barrierSfii 
strong intramolecular correlations/^ band-gap inducements^ 
and so on. One of early experiments by Li et al suggested that 
as the NDR region was approaching, the suppressed value of 
the Fano factor would increase,^ Further exploration by lan- 
naccone et al^ showed that in the NDR region, the shot noise 
would go through a transition from a sub-Poissonian value 



to a super-Poissonian value in a nonlinear fashion. Never- 
theless, a NDR was not a sufficient condition to generate the 
enhancement. As Song et al showed that there was no noise 
escalation in the super-lattice tunnel diode even though its I-V 
curve also exhibited a NDR mg\onXL This led to a conclu- 
sion that charge accumulation, which was related to the in- 
ternal Coulomb potential, was ultimately responsible for the 
super-Poissonian shot noise. Given the good agreement be- 
tween numerical calculations from semi-classical theory and 
those experiments, the Coulomb interaction was thought to be 
the reason for the noise enhancement. Since quantum effect 
dominates the transport behavior in mesoscopic systems, a 
quantum theory of shot noise capable describing the enhance- 
ment in the NDR region is clearly needed. In 1999, Blanter 
and Buttiker have studied the shot noise of resonant tunneling 
quantum well theoretically using scattering matrix method 
In the nonlinear regime. Coulomb interaction (Hartree level) 
leads to hysteretic behavior in I-V curve. By including the 
fluctuation of Coulomb interaction, they identified an impor- 
tant energy scale, termed interaction energy, in the Fano fac- 
tor. They found that in the NDR region where the interaction 
energy is very large, a super-poissonian behavior occurs due 
to the fluctuation of Coulomb interaction. 

Understanding electronic transport properties of atomic- 
wire based structures is very important from the scientific 
viewpoint and due to its potential applications in molec- 
ular electronics. For example, combining the Lippmann- 
Schwinger equation and density functional theory (DFT), a 
NDR in the tunneling regime of atomic carbon wires was 
predicted by Langji^ The shot noise of silicon atomic wires 
has also been studied using the same approachaS. In this 
paper, we develop a general theory for dc shot noise by in- 
cluding the fluctuation of Coulomb interaction. Our theory is 
based on non-equilibrium Green's function (NEGF) method 
which can be coupled with DFT to study transport properties 
of nano-devices from first principles. As an application of 
our theory, we investigate the shot noise of an atomic carbon 
wire structure with four carbon atoms in the scattering region 
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(Al — C4 — Al). Its I-V characteristic and transport proper- 
ties have been well understoodi^ It was found that a band gap 
induced NDR occurs at high bias due to a shift of conduction 
channels in the central region. We have used the traditional 
formula2ai2^ 

Sa^p = (l/2)[< At(i)A/^(t') > + < Aip{t')AL{t) >] 



dE{[fM - fa) + Ml - M]Tr[f] 



(1) 



to calculate the shot noise for Al — C4 — Al structure. Our 
results show that shot noise is sub-poissonian. When the fluc- 
tuation of Coulomb interaction is included, large shot noise 
was found in NDR region showing super-poissonian behav- 
ior. 

Our paper is organized as follows. In section II, we de- 
rive a general theory for dc shot noise when the fluctuation of 
Coulomb interaction is included in the first order The detailed 
derivation is given in Appendix. In section III, we describe 
some technical details and show the numerical results in the 
atomic carbon chain system along with an analysis and dis- 
cussion of the result. Finally, the summary is given in section 
IV. 



n. THEORETICAL FORMALISM 

In this section, a NEGF theory is developed to calculate dc 
shot noise in the regime of NDR, which involves the Coulomb 
interaction between electrons. A key ingredient of the new 
theory is that to account for large shot noise in the NDR re- 
gion both self-consistent Coulomb potential and its fluctuation 
have to be considered. 



A. General Expression 

We start from a quantum coherent two-lead conductor de- 
fined by the Hamiltonian 

ka n 

+ Y.\-*kanCLdn+C.C.] (2) 

kan 

where Cl^{Cka), dl^{dn) are the creation (annihilation) op- 
erators of electrons in leads and the scattering region, respec- 
tively. The first term describes leads that dc voltages are ap- 

phed on, and e^a = ^fl + l^a which e^*^^ is energy levels in 
the lead a and Va stands for an external voltage. The sec- 
ond term is for the isolated central region, where the self- 
consistent internal Coulomb potential under the Hartree ap- 
proximation is defined as 



Un 



where Vnm is a matrix element of the Coulomb potential. In 
the real space V{x, x') = l/\x — x'\, q is the electron charge. 
The last term corresponds to a coupling between the central 
region and leads described by a coupling constant tkan ■ 
The current operator of the lead a is defined as (h = 1): 



(4) 



where Na = J2k ^ka^ka is the number operator for elec- 
trons in the lead a. 

From the Heisenberg equation of motion. 



"df 



we have 



dNo, 
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'■^[tkanCladn] + h.C 



(5) 



(6) 



kn 



(3) 



where H is the system Hamiltonian with constant Coulomb 
potential and h.c. denotes the Hermitian conjugate. 
Hence the current operator becomes 

/a„(t) = -iqY,[tkanCijt)dnit)] + h.C. (7) 
kn 

On the mean field level, the current is a functional of 
Coulomb interaction, i.e., / = /[([/)]. Here we have treated 
the operator of Coulomb interaction [/ as a C number mean- 
ing that the fluctuation of Coulomb interaction is assumed not 
important. However, this is not always true. For instance, in 
order to reflect the Coulomb interaction between electrons in 
the NDR region, we have to consider the fluctuation of the 
Coulomb potential. Fig.l shows the physical picture of the 
NDR. For simplicity, we assume that the scattering region has 
one resonant level Eq with a width characterizing the lifetime 
of the resonant level. We also assume that there is an effec- 
tive band bottom for the lead which is crucial for the phe- 
nomenon of NDR. As shown in Fig.l, when the bias voltage 
is increased the current increases because the resonant level is 
brought down by the external bias. As the bias is increased 
further such that the resonant level falls below the band bot- 
tom of the lead the current starts to decrease giving rise to the 
NDR. The above physical picture is static where the Coulomb 
potential is included on the mean field level and the correla- 
tion effect of Coulomb interaction has been neglected. For the 
current correlation in the NDR region, the correlation effect 
of Coulomb interaction has to be considered. In this picture, 
when the resonant energy level Eq is about to fall below the 
band bottom of the lead, the internal potential of the scattering 
region due to the Coulomb interaction of injected electron will 
push it up, leading to a positive correlation between incoming 
electron flows. This positive correlation is a dynamic process 
and can not be described by a Hartree field. In another word, 
the fluctuation of Coulomb interaction has to be considered for 
the positive correlation in the NDR region. As demonstrated 
by LaradeSi and confirmed by our calculation, for the atomic 
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FIG. 1 : Graphical band profile for nanoscale devices in the presence 
of bias. At the zero bias, the resonant level is above the Fermi level 
(a). As the bias increases, the Fermi level of left lead increases and 
is above the resonant level (b) and (c), giving rise to a large increase 
of current. As the bias increases further, the resonant level is below 
the band bottom of the left lead (d) leading to a NDR. 



wire with even number carbon atoms like C4 and Cg, there is 
an effective band bottom responsible for the NDR. Therefore, 
the fluctuation of Coulomb potential U should be important. 
For odd number wires such as C5 and C7, however, there is no 
apparent NDR effect. Hence there no effective band bottom 
and the fluctuation of U can be neglected. 

To treat the fluctuation of Coulomb interaction, we follow 
the idea of Reffisl and expand the cuiTent in terms of Coulomb 
potential operator about its equilibrium value up to linear or- 
der After the expansion, the total current in the real space 
could be expressed as 



6Ia 

su,, 



(8) 



where /„ =< ia„ >, Ui =< Ui >, and SUi{t) = Ui{t) - J7,. 
We have also introduced a quantity Xai = Sla/SUi. 

It is easy to see that this new current operator gives the same 
current but different shot noise. With the new current operator, 
the current correlation is obtained up to the linear order in SU, 



sl^> AUt)Aip{t') >=< {L{t)ip{t') > -i^ip 



(9) 



where 3^/3 = (1/2) [S^^^+S^^^]. The first two terms in Eq.® 

corresponds to the current correlation sj^^^ in the absence of 
Coulomb potential fluctuation. They have been calculated be- 
fore and is given by Eq.([T]i. The last term in Eq.(|9]l is defined 



^(1) 



i 

-{K^uap + Xp^ujo.)] (10) 



Using the NEGF method, we have derived the expression of 

(1) 



(see the appendix). Finally, the shot noise in the presence 

(11) 



of Coulomb potential fluctuation is written as 
where 



(12) 



with 

Xai 

and 



-j- [ dEY^ifa - MiGT^G^TaG''),, + h.c. 

(13) 

fi£;[G"^r„G''S>G"'7a - G"'S<G°r„G°(/„ - 1) 

+ G''S]<G°r„G''S]>G°] (14) 

Eqs. (fTT]) - (fT4l) is the central result of this paper. To calculate 
the shot noise, one has to solve the Green's function together 
with the Poisson equation. 



G'' 



1 



E-H-U-V,c 



and 



V^[/(x) = Aniq 



dE 
2^ 



G<{E,x,x) 



(15) 



(16) 



where V^c is the potential due to the exchange and correlation 
effect in the first principles calculation. 

For a two-terminal device at zero temperature, we set a = 
L,f3 ~ R and vu > vl = so that — 1) = and 

fnifL — 1) = at zero temperature. We then obtain 



; J dE[GTLG''T.>G''fL 



- G''Y.<G''TLG°-{fL - 1) + G''S<G"rLG'^I]>G''] 

dEiCTLC'TRG" + iGTLGTLGT^G'') 

(17) 

Similarly, we have 



Ep-qvR 



dEiGTRGTLG'' + iGTRG'TRGTLG") 

(18) 
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Lead Scattering Region Lead 



FIG. 2: The schematic plot of Al — C4 — Al system. The atomic 
wire with four carbon atoms is linked by two semi-infinite Al elec- 
trodes. The (100) direction Al electrodes extend to ±00 where elec- 
tron reservoirs are located. 

In addition, from Eq.(fT3]|. we have 

Al = - ^ / " dE{G'TRG''TLG') + h.c. 

>^R^^ f ' dSlG'^FiG'T^G'-) + h.c. (19) 

ZTT JEp-qvn. 

With Eqs.dTTli. ( fTSl l. and ( fT9] l. the two-probe shot noise can be 
calculated. 

III. NUMERICAL RESULTS 

In this section we use state-of-the-art first-principle quan- 
tum transport package MATDCAL to investigate the general 
transport properties of atomic carbon-chain systems coupled 
with Al leads. In the package DFT is carried out within 
the formalism of the Keldysh nonequilibrium Green's func- 
tion. Numerically, the effective Kohn-Sham (KS) equations is 
solved by a linear combination of the atomic orbitals (LCAO) 
basis set. We define the atomic core by a nonlocal norm con- 
serving pseudopotential and treat the exchange-correlation at 
the LDA level. DFT determines the atomic structure and the 
system Hamiltonian while NEGF contributes to the nonequi- 
librium transport properties. Under an external bias the trans- 
port boundary conditions are treated by the real space numer- 
ical techniques. For further references, the theoretical back- 
ground and practical execution of this formalism can be found 
in Refi24i A numerical error tolerance is set to be 10^"' to 
confirm self-consistency. 

Generally speaking, we have carried out our calculation on 
the atomic chain structure with four carbon atoms Al — C4 — 
Al^ The carbon chain lies in the central simulation box in 
contact with electron reservoirs through two semi-infinite Al 
electrodes. The schematic structure is shown in Fig. 2 where 
there are 18 Al atoms in the unit cell of the semi-infinite elec- 
trodes with a cross section along (100) direction. The con- 
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Voltage 

no. 3: The I-V curve of Al - d - Al structure. A NDR region 
begins to show up around 0.65 V. Inset: the transmission coefficient 
for this system at zero bias. 

tact distance between the Al electrode and the carbon chain 
is fixed at 0.378 a.u. while the distance between the nearby 
carbon atoms is equal to 2.50 a.u. In our calculation, we have 
set temperature to be zero. 

Technically, the correction term in Eq. (fT2l) has to be solved 
in real space due to the Coulomb like interaction involving 
Vij which reads V{x,x') = 1/|.t — x'\ in real space. Since 
quantities A and S play the role of charge, we can define the 
potential induced by S, 

^ax = '^V{x,x)lTSl{^a{x' ,x')) (20) 

x' 

Since V^Y {x,x') = —4:1t5{x — x'), ila satisfies the Poisson 
like equation, 

y^n^ix) = ~ATTlm{E^{x,x)) (21) 

We solve this equation for the leads to obtain the boundary 
condition for the scattering region. It turns out ila{x) in the 
lead is very small so that the boundary condition of fl^ {x) can 
be safely set to zero. Once fla (x) is obtained the correction 
term can be easily calculated from Eq.(fT2]i 

g2 f 

^aP = / [>'l3{x)^aix) + Xa{x)V,i3{x)]dx (22) 

The I-V characteristics is shown in Fig. 3, where the inset 
plots the transmission coefficient T versus the energy E at zero 
bias. The shot noise and the corresponding Fano factor are 
shown in Fig. 4. In Fig. 3 and 4, the Coulomb interaction is 
included on the Hartree level and the Coulomb potential fluc- 
tuation is neglected. Following observations are in order. (1) 
The I-V curve is similar to the result obtained by Larade et 
ali^ (2) At zero bias the resonant energy is higher than Fermi 
energy of the system (chosen to be zero). As we apply a volt- 
age to the right electrode, the resonant energy level should 
drop and move closer to Fermi energy of the left electrode. 
With the increasing of voltage, the effective band bottom of 
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the emitter and the main resonant level would be aligned, and 
this gives the maximum current around 0.65 V. When the volt- 
age increases further, a significant decrease of current occurs. 
(3) The shot noise in the absence of Coulomb potential fluc- 
tuation has a similar behavior as that of I-V except that the 
maximum is at 0.7V instead of 0.65V. (4) The Fano factor is 
nearly a constant of order of 0.6 in positive differential resis- 
tance (PDR) region when bias is smaller than 0.6V. It starts to 
increase sharply upon entering the NDR region and the Fano 
factor eventually shoots up to 0.8. We conclude that although 
the Fano factor calculated on the Hartree level shows enhance- 
ment in the NDR region, it is still sub-Poissonian which does 
not agree with experimental resulti^ 

In general, when electrons tunnel through the left barrier 
to occupy empty energy levels, Pauli principle inhibits other 
tunneling electrons to reach the same energy level but higher 
ones. As a consequence Pauli exclusion principle gives the 
negative effect for current correlation. Coulomb interaction, 
however, can give either positive or negative correlation effect 
to shot noise. This can be understood as follows. It is known 
that the maximum current corresponds to the situation that the 
energy of incoming electron is in line with the resonant level. 
Hence in both PDR and NDR regions, most of electrons are 
off resonance. Due to the Coulomb interaction the incoming 
electron can push up the resonant level so that the next elec- 
tron will be further away from the resonance in the PDR re- 
gion or closer to the resonance in the NDR region, giving rise 
negative or positive correlations. Our numerical results indeed 
confirm this physical picture. In Fig. 5 and 6 we present the re- 
sult of shot noise and Fano factor in the presence of fluctuation 
of Coulomb interaction. We have also included the shot noise 
and Fano factor in the absence of Coulomb potential fluctua- 
tion for comparison. In Fig. 5, we see that the correction term 
solved via the Poisson like equation is very small at low bias 
when V < 0.5 and becomes negative until around l.OV where 
shot noise increases sharply. In Fig. 6, we see that a large Fano 
factor great than 3 occurs near V — 1.0. This result is in 
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Voltage (V) 

FIG. 6: The Fano factor plotted as a function of the apphed voltage 
for Al — Ci — Al system at K. Circle: the Fano factor is calculated 
by the Eq.(|TJ for the same system. 

qualitatively agreement with others' work i^-'^'^^ 

IV. CONCLUSION 

The traditional formalism can only describe the suppres- 
sion of shot noise, which corresponds to the PDR region. In 
order to treat enhancement of shot noise correctly in the NDR 
region, we have to include the fluctuation of Coulomb interac- 
tion. In this paper we have developed such a general dc theory 
for calculating the shot-noise in the NDR region. The theoret- 
ical framework is based on the combination of NEGF-DFT 
formalism with the self-consistent Coulomb potential and its 
fluctuation included. Our theory (Eq.(fTTTi) can also be applied 
to mesoscopic conductors. We have applied our theory to 
molecular devices. Specifically, we have calculated the shot 
noise of AI — C4 — Al structure which is an ideal system since 
its I-V curve exhibits a NDR region. We found a large Fano 
factor in the NDR region exhibiting super-Poissonian behav- 



ior. 



which is equivalent to Eq.fTsTl. 
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V. APPENDIX 

In this appendix, we will derive the expression of Xa^p and 
Aq,^ using the theory of NEGF. 

Expression of Aq ,9 

Using the NEGF, the current is given by 

/a = / di?^[rr(r„GT/3G")](/„ - fp) (23) 

where f^.p are Fermi distribution functions in corresponding 
leads, G""'° are the respective retarded and advanced Green 
functions of the scattering region and Fq ,g are the linewidth 
functions related to the coupling of leads and the scattering 
region. 

To obtain Aq, we first calculate the following quantity, 

d 



(24) 



The quantity duG^ can be calculated from the Dyson 



equations^ 



G'' = G'o + GlUG'' 



(25) 



where Gq is the retarded Green's function in the absence of 
Coulomb interaction. Taking the derivative with respect to 
Ui,we have 

8 8 8C^ 

—G- = —{GlUG')^GlD,G- + GlU— (26) 

where Di is a diagonal matrix with the matrix element 
{Di)jk = SjiSki, i.e., there is only one nonzero matrix ele- 
ment. From Eq.([26]l, we find 



_d_ 



G'' = 



1 



1 - GIU 



-GlDiG'' = CD^C 



(27) 



where we have used the following Dyson equation again, 

1 



G'' 



1 - GIU 



Substituting Eq.(l27]i into Eq.(l24li. we can get 

^o.^iE) = E rr(r„G''AG'^r;3G") + h.c. 



2tt 



E Tr(G'T0G°r„G''A) + h.c. 



-T-E(G''r^G"r„G^)., + ;i.c. (28) 



Expression of A 



a/3 



For simplicity, we only deal with the first term of Aap ex- 
plicitly, i.e. < Iaa{t)Ui{t') >. The second term can be cal- 
culated similarly. Using the current and Coulomb potential 
operators in Eqs.(|7]i and (O, we obtain 

< —iq 

kn m 
kmn 

-Cn<4WC'fec.(t)4(i')4.(i') >] 

= E < ^L(*)rfm(i') >< d„{t)dlit') > 



C„ < di{t)d^{t') X Gfc„(t)4(i') >]+IM 



(29) 



In terms of Green's functionsai, 

Gnm{t',t) = -l<dn{t')dl{t) > 
G<Jt,t')^i<di{t')d^{t) > 
GLmit',t)--^<CUt')dl{t)> 



(30) 
(31) 
(32) 
(33) 



we obtain. 



= -iq^ E ^^^n[tkanG'^ i.^{t,t')G^^^{t',t) 



knm 



(34) 



Applying the Langreth theorem of analytic continuation^^ and 
suppressing time indices, we have 

Gm,ka = ^'^{Gmltkal9ka + G^it'l^ig^^) (35) 
Gka,m = ^i9kJkalGi„, + oUkalGf^J (36) 



where g^^'^'"^ are the corresponding Green's functions in the 
lead a. 

For dc transport, the Green functions depend only on t' — t. 
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After the Fourier transform from time to energy, it becomes 

knml 

~^ G ml^ kal9 ka)G nm ~ ^kan^mnidka^kalGi^ 

+ 9kJkalGf^n)] 
■ 2 r 



we find 



+ S]q,G^ )]mm 

^""J 



where ''''' are the corresponding functions of self-energy 
due to the lead a. Using the Keldysh equation and properties 
of Green's functions^, 

G< = G''I]<G'' (38) 

Y.l{E)-Y;'^{E)=iT^{E-qv^) (39) 

S<(S) = iT^{E - qv^)UE) (40) 

E>(£)=zr„(^-gt.„)(/„(S)-l) (41) 



S„ = I j dE[G'V^G^Y>G''f^ - G^Y.<G''T^G%f^ - 1) 
+ G''S<G°r„G''S>G°] 



(42) 



Substituting this equation into Eqs.dJTl) and ( fTOl i. we finally 
arrive at 



= E(V^^^^-2a. - Aa.F.,S^,) (43) 



(37) 

where we have used the fact that < Iao{t)Ui{t') >=< 



U^{t)I^,{t') >T. 
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